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1 Introduction 


We present a new algorithm 

n 

A n(f) = '^Ta i f(x i ) 
i=l 

for the approximation of integrals 

W) = / Six) Ax. 

J[ o,i] d 

Fred Hickernell wrote a paper “My dream quadrature rule” where he proposed 
five criteria that an ideal or “dream” quadrature formula should satisfy. We also 
present a list of five (similar, but different) properties of our “dream algorithm”: 

(pi) The algorithm A n should be an unbiased randomized algorithm, i.e., 

E (MS)) = W) 

for all integrable functions. Of course this means that the weights a* G M. and 
the points Xi G [0, l] rf are random variables. It is beneficial to have positive 
weights at > 0 for all i. 

(P2) The randomized error 

E(K(/)-/ d (/)|) 

of A n should be small and/or optimal in the sense of order of convergence 
for “many” different classes of functions. In particular, we would like to have 

E(\A n (f) - I d (f) |) < c r4 n~ r ~ l/2 (logn) (d-1)/2 ||/||^.-- ([ o i i ] d ) (1) 

for all r G N, as well as for all sGN with s > d/2 

E(|A,(/) - IiU)\) < II/II».([0,1K). (2) 

(P3) The worst case error 

supK(/)-/ d (/)| 

LO 

among the realizations Ad/ of A n should be small and/or optimal in the sense 
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of order of convergence for “many” different classes of functions, in particular 


SU P(\ A nU) - W)\) < c r,d n T ( lo S n) [d 1)/2 ||/||iJi-,mix ([0)1] d ) (3) 

UJ 

for all r 6 N, as well as for all s G N with s > d/2 

sM\K(f) - W)\ ) < c s , d n~ s/d ||/||h.([ 0 ,!]<«)■ (4) 

(JJ 

(P4) The algorithm should have good tractability properties in the sense of the 
theory of “tractability of multivariate problems”, see [5]. 

(P5) The algorithm should be easy to implement. 

In this paper we concentrate on properties (PI), (P2) and (P3) and hence 
we are not specific on (P4) and (P5) and leave them for further research. In 
particular, we do not discuss tractability and all constants c > 0 may depend on 
the dimension d and the smoothness r or s. A few remarks are in order. 

1. The simplest Monte Carlo method certainly satisfies (PI). Therefore it is 
easy to run the algorithm a few times and to do an (a posteriori) error 
analysis. This is a great advantage of an unbiased algorithm. Of course the 
low rate n~ 1 ' 2 (even for very smooth integrands) is a big disadvantage of 
the simplest Monte Carlo method. Randomized algorithms with a higher 
rate of convergence are known and often they are unbiased; usually they are 
designed for a specific class of functions. 

2. We do not know of any algorithm in the literature that satisfies (P2), even 
in the univariate case d = 1. The upper bound ([[]) seems to be new. The 
main term n~ r ~ 1//2 is of course optimal. 

The bounds (2)-(4) are known and it is also known that they are optimal. 
The bound (3) is from Frolov, see [HEIE]. The bounds (2) and (4) are from 
Bakhvalov and can be found in [4|. 

3. Many known algorithms (such as the Gaussian quadrature formulas) satisfy 
(P3) in the univariate case. It is also known that (modifications of) the 
Frolov algorithm satisfy (P3) for arbitrary d. Hence the Frolov algorithm 
(or some modifications of it) is “universal” in the worst case setting, see also 
the recent paper nnj. Since it is a deterministic algorithm it certainly cannot 


3 


satisfy (PI) or (P2). The problem with any deterministic algorithm A n is 
that a computation of A n (f) does not come together with an error bound 
since usually the norm of / is not known. 

4. We did not discuss the property “extensible” in the list of Hickernell. We 
believe that this is another nice property but not as important as the other 
properties since it can decrease the total computing time only slightly. 

In this paper we present an algorithm M a ^s with positive weights that satisfies 
(pi) and (P2) and (P3), see Section 0 In particular we prove the existence of 
A n such that ([]]) holds. 


2 Some Notation 


For r, d G N the tensor product Sobolev space H r,mix (R d ) is defined as the space 
H r ' mix (R d ) = {/ G L 2 (R d ) | D a f G L 2 (R d ) for every a G {0,.. .,r} d J 
of real valued functions, equipped with the scalar product 


(Ad) 


(JJd) 


E < B "/. B ”9W) 

a£{0,...,r} d 


and hence with the norm 


fjr, mix 


1/2 


E ii D "/n 


L 2 (K d ) 




It is known that H r,mix (R d ) is a Hilbert space and its elements can be taken to be 
continuous functions. In this paper, the Fourier transform is the unique continuous 
linear map : : L 2 (R d ) — > L 2 (R d ) with 

/(</)= / f(x)e- M ^Ax 
J R d 

for integrable / and y G R d . The space H r ’ mix (R d ) contains exactly those functions 
/ G L 2 (M d ) with / • }\J 2 G L 2 (M d ) for the Fourier transform f of f and the weight 
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function 


d r 




", h r (x)= y = nY\ 2nx i 

ae { 0,...,r> d J=1 j=l k=0 


2k 


In terms of its Fourier transform, the norm of / G H 


r,mix/Tn)e£\ 


is given by 


2 

fjr, mix 


f(x) ■ h r (x) dx. 


Analogously, for s, d G N the isotropic Sobolev space H s (R d ) is the space 
H s (R d ) = {/ G L 2 (M d ) I D a f G L 2 (R d ) for every a G with ||a||i < s} 
of real valued functions, equipped with the scalar product 


if, 9) 


H‘ 


E (D“f,D°g) 


l 2 ( m d ) 


a l<s 


and hence with the norm 


1/2 


H‘ 


E ii^/ii 


L 2 (R d ) 


\\a i<s 


This also defines a Hilbert space. In the following, let s > d/2. Then H s (R d ) 
also consists of continuous functions, exactly those functions / G L' 2 (R d ) with 
/ • v l J 2 G L 2 {R d ) for the Fourier transform / of / and the weight function 


v. : -»■ 


^ x = 


e ni^-i^-a+iNih 


IM|i<s i=l 


In terms of its Fourier transform, the norm of / G is given by 


/(x) -i; s (x)dx. 


Furthermore let C' c (M a! ) be the set of all continuous real valued functions with 
compact support in R d . 

We will first present an unbiased Monte Carlo method for integration on (7 c (M d ) 
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in Section [4j We will examine its error for the subspaces H r,mix (M. d ) and H s (R d ) of 
functions in H r,mix (R d ) or H s (R d ) with compact support. This includes an error 
bound for the classes H r,mix ([0, l] d ) and H s ([ 0, l] d ) of all functions in H r,mix (R d ) or 
H s (M. d ) with support in the unit cube [0, l] d . These spaces can also be considered 
as subspaces of the Hilbert space 

TT’ mix ([0, l} d ) = {/ e L 2 ([0, l] d ) I D a f e L 2 {[ 0, l] d ) for every a e {0,... ,r} d } , 
equipped with the scalar product 

(/) #)m->mixQ 0 ,l]d) = (D f,D 5 , )i2([ 01 ]d) , 

ae{0,...,r} d 


or the Hilbert space 

# s ([0, l] d ) = {fe L 2 ([ 0, l] d ) I D a f e L 2 ([ 0, l} d ) for a G No with \\a\h < s} , 
with the scalar product 

(/)<7)h«([ o,i] rf ) = (D a f,D a g) L 2([ o,ip), 

||a||i<s 

respectively. It turns out that this method for // r,mix ([0, l] d ) and FF([0, l] d ) can 
be transformed to the full spaces hP' ,mix ([0, l] d ) and H s ([ 0, l] d ) without loosing its 
good properties. 

3 The Basic Quadrature Rule Qs iV 

Let S e be any invertible matrix and v any vector in At the basis of the 
Monte Carlo methods to be presented is the deterministic and linear quadrature 
rule Qs,vi defined by 

QsM) = isy E 1 ( s " (m +” 

me z d 

for any admissible input function / : W l —> M. This includes all functions / 
with compact support. For such functions the sum is actually a finite sum. More 
precisely, Qs, v uses the nodes S~ T (m + v), where m e 7L d is a lattice point in the 
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compact set (S' 7 (supp /) — w) of Lebesgue measure (det(S') • (supp /)). This 
volume is the approximate number of nodes of Qs, v - 

In particular, the number of nodes of Q a s,v for a > 1 is of order a d . The 
following simple lemma gives an exact upper bound, see [6] for other bounds. 


Lemma 1. Suppose f : M. d —> M. is supported in an axis-parallel cube of edge length 
l > 0. For any invertible matrix S E v E W l and a > 1 the quadrature rule 

QaS,v uses at most (Z - ||-S'||i + l) rf - function values of f. 


Proof. By assumption, / has compact support in | • [— 1 , l^ + xo for some Xq E 
The number of function values is bounded by the size of 

M = |m E 7L d | (aS)~ T (m + v) E ^ • [—1, l] rf + x 0 
= im E Jf' \ m + (y — aS T x o) E — ■ S' 7 !—1, l] d ^ . 


Since ||S' T x|| 00 < ||S' 7 || 00 = ||S||i for x E [-1, l] d , 
MC l m E Z d \ m + (y — aS T x 0 ) E 


fl|S||i,fl|S||i 


and \M\ < (a/||S'||i + l) d . With 1 < a we get the estimate of Lemma [U □ 

The error of this algorithm for integration on C c (R d ) can be expressed in terms 
of the Fourier transform. 

Lemma 2. For any invertible matrix S E R dxd , v E and f E C c (M, d ) 


\Qs.M-W)\< E 


f(Sm ) 


m£ Z d \{0} 

Proof. The function g = f o S~ T (■ +v ) is continuous with compact support. Hence, 
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the Poisson summation formula and an affine linear substitution x = S T y — v yield 


QsM) — 




m = 


i detS i„s< l detS l 


s( m ) 


\ieTs]^J f{S ' T{X + V)) 

me* Dd 


. e ~2m(x,m) ^ 


f f(y) 

nez d md 


[1 —2iri{S T y—v,m) 


d y 


f(Sm) ■ e 2 ™ {v ' m) , 


if the latter series converges absolutely, see [31 pp. 356]. If not, the stated inequality 
is obvious. This proves the statement, since /<*(/) — f(S ■ 0) • e 27ri ^ ,0 >. □ 


4 The Method M a ^ for Integration on H r,TDJX (R d ) 
and H s (M d ) 


It is known how to choose S in Qs, v to get a good deterministic quadrature rule 
on iL r,mix ([0, l] d ). Let the matrix B G M, dxd satisfy the following three conditions: 

(a) B is invertible, 


(b) 


d 


n ( Bm )j 

3 =1 


> 1, for any m G 7L d \ {0}, 


(c) For any x, y G W l the box [x, y\ with volume V 
most V + 1 lattice points Bm, m G 


d 

[ \xj — yj | contains at 

3 = 1 


where [x, y] = {z G | Zj is inbetween of x 3 and yj for j = 1,..., dj. Such a ma¬ 
trix shall be called a Frolov matrix. Property (b) says that for a > 0 every point 
of the lattice aBZ d but zero lies in the set D a of all x G with n-=i \xj I > a d . 







This graphic shows the lat¬ 
tice aBTJ 1 for d = 2, a = 3 
and the Frolov matrix 



B 


1 2-V2\ 

1 2 + y/2) 


Except zero, every lattice 
point lies inside D ?1 . 


It is known that one can construct such a matrix B in the following way. Let 
p E 7 j[x\ be a polynomial of degree d with leading coefficient 1 which is irreducible 
over Q and has d different real roots Ci, ■ ■ ■, Ci- Then the matrix 


B 



d 


i,j =1 


has the desired properties, as shown in |3 p. 364] and |9]. In arbitrary dimension 
d we can choose p(x) = (x — l)(x — 3) • ... • [x — 2d + 1) — 1, see |T] or [9], 
but there are many other possible choices. For example, if d is a power of two, 
we can set p{x) = 2 cos (d ■ arccos(x/2)) = 2Ta(x/2), where T ( i is the Chebyshev 
polynomial of degree d, see [7] p. 365]. Then the roots of p are explicitly given by 
C j = 2 COS (^y-vr) for j = 1,..., d. 

K. K. Frolov has already seen in 1976 that the algorithm Q a B,o for a > 1 is 
optimal on LT r,mix ([0, l] d ) in the sense of order of convergence. It satisfies 

\QaB,o(f) ~ T d(f )I < ca~ rd ■ (logo) V • ||/|| Hr , mix([0il]<J) 

for a constant c > 0 and any a > 2 and / G H r,mix ([ 0, l] d ). We hence call it Frolov 
quadrature formula. See also [I] or [9] for a proof. In fact, the same error bound 
holds for QaB,v for any v G M d . 

We define a randomized version of this quadrature rule by introducing two 
independent random vectors v and u. With the random shift parameter v G [0, l] d 
the algorithm gets unbiased. The random dilation parameter u G [1,2 1 / rf ] rf will 
ensure the general error bound of Theorem 1. Both effects are independent of 
each other: The random shift is not needed for the error bound and the dilation 
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is not needed for the unbiasedness. 

Algorithm. For a Frolov matrix B G M, dxd and any a > 0 the randomized Frolov 
quadrature formula M aB is the method Q a uB,v from Section [3] with independent 
random vectors u and v, uniformly distributed in [1, 2 1 ^ d ] d and [0, l] d respectively 
and u = diag(rti,..., ufi). 

Lemma 3. The method M a ^ B is well-defined and unbiased on L 1 (R d ). 

Proof. We realize that for / G L l (W l ) 

\QauB, v (f)\ Y” , /-pi / |/ ((auB)- T [m + x)) \ dx 

^ d \detauB\J md ' 

= E u [ \f(y)\ dy 

m&d J{auB)-r m +{auB)-T[0,l]d 

= Eu / |/(l/)I dy = I \f(y)\ dy < oo. 

J JR d 

We can thus apply Fubini’s theorem and get 

E(Ma, B (f)) =^uK(QauB,v(f)) 

= E « u i. 1 — r>| l f {{auB)- T (m + x)) dx 
Idet auB\ J [0A]d 

= E u V f f{y) dy 

m£ Z d BauB)~ T m+(auB)~ T [ 0,l] d 

= E « [ f(.y)dy = W)- 
J R d 

In particular, M atB (f) is almost surely finite. □ 

According to Lemma[T]the method M a B uses no more than 2-(l ■ ||-B||i + 1 ) d -a d 
function values of a function / supported in a cube of edge length l. Later we will 
show that M a ^ B satisfies 

E|- «/)| < ca~ rd ~ d ^ • (log a)^ ■ ||/|| H „„,. (RJ) 

for a constant c > 0 and any a > 2 1 / d and / G H r,mix (M. d ). But first we analyze 
M a)B on the larger set C c (R d ). 
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Error Bound for C c (W l ) 

We prove a main result of this paper. Again, D a is the set of all x G with 
UU \ x o\ — The me fh°d M a B satishes a general error bound on 

Theorem 1. Let B G M dxd 6e a Frolov matrix. Then there is a constant c > 0 
such that for every a > 0 and f G 


E| M a>B (f) 


W)\ < ca~ d 



dx. 


Proof. Let v G M. d be arbitrary, but hxed. Thanks to Lemma [2] and the monotone 
convergence theorem we have 


Eu I QausAf) - W)\ < E u [ f^uBm) 

, mez d \{0} 


E. u f(auBm) 

mGZ d \{ 0} 


Since each auBm is uniformly distributed in the box [aBm, 2 1 / d aBm] with volume 
(2 1 / d -l) d - uU a(Bm)j\, this series equals 


_ - _ V" 

( 2Vd - !)' me ^ { o } 
1 


f( x ) 


[aBm.,2 1 / d aBm] 


n;=i H Bm )o 


dx 


< 


( 21/d - l f m Jt<0) 


fi x ) 


[aBm,2 1 / d aSm] 

f{x 


n,,i 2 - i/d 


dx 


Xj 


(2Vd - l) c 
2 


n u to 

fi x 


(2 1 / d — l) d V UU \ x i 


{ m G h d \ {0} | x G [aBm, 2 lB aBm]} | dx 

, U/d \ I 


Bm G 

‘ X X' 

) 


i2 l / d a J a\ 

/ 


dx. 


Thanks to the properties of the Frolov matrix B, if fXj=\ \ x j\ < ad i d ie latter se t is 
empty and otherwise contains no more than fljLi |^f| + 1 < 2 a~ d rij=i I x j\ points. 
Thus, we arrive at 


Vu\QauB,v(f) ~ W)\ < 


„ —d 


(21/d _ j Da 


fi x ) 


dx. 
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By Fubini’s theorem, we have 


E | M a>B (f) - I d (f)\ = E„E U | QausM) - W)\ < 


-d 


(2 1 /* - l) c 


>D a 


fix) 


and the theorem is proven. 


dx 


□ 


R d ) result in decay properties 


dx. Hence, the general 


Additional differentiability properties of / G G c 
of /. This leads to estimates of the integral f D f(x) 
upper bound for the error of M a ^{f) in Theorem Q] adjusts to the differentiability 
of /. Two such examples are functions from H r,mix (R d ) and H s (M. d ). 


Error Bounds for H r ' mix (R d ) 

If / G H r,mix (R d ) C C' c (M ri ), the following lemma holds. 


Lemma 4. For any Frolov matrix B G R dxd and r G N there is some c > 0 such 
that for each a > 2 1 / d and f G H r,mix (R d ) 


'D a 


/(x) d x < ca rd+d / 2 (logo) 2 


fjr, mix 


Proof. Applying I lolder’s inequality and a linear substitution x = aBy to the 
above integral, we get 


'D a 


< 


fix) 


2 

fjr, mix 


dx ) = 


'D a 


f(x) hr(x) 1 ^ 2 ■ h r (x) ^ 2 dx 


h r (x) 1 dx = 


H r 


'Da 


'G 


h r (aBy) 1 d^/ ) - or |det B \, 


where G = B~ 1 D\ is the set of all y G with f|j=i \ {By)j\ > 1. It it thus 
sufficient to prove that the integral J G h r (aBy)^ 1 dy is bounded by a constant 
multiple of a~ 2rd ■ (loga) d_1 . 

Consider the auxiliary set N((3) = {x G R d | < \xj\ < 2^, 1 < j < d} for 

P G Nq. Let \P\ = EUiPi and G « = G n {v G | aBy G N(P)}. Since R d is 
the disjoint union of all N(P), G is the disjoint union of all Gf over P G Nq. For 


y G Gf we have both Mill a(By) 


> a d , since y G G, and 


If ,«(%), < 2^1, 


since 


aBy G N(P). This implies Gf = 0 for \P\ < dlog 2 a, since then 2^1 < a d . 
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Let y e Gf and \(3\ > dlog 2 a. Then 

hriflBy) > JJ (l + L2^- 1 J 2r ) > Jj2 2r( ^'- 1) = 2 2r( l /3| - d) 

J=1 j =1 

and hence h r (a5y) _1 < 2 2 b ( H£l). On the other hand 

A d (Gf) < A, ((ai?)- 1 ^^)) =a~ d -\ det B\~ l ■ \ d (N(P)) 

d 

= a~ d ■ | det 5| _1 2 d Yl ( 2 ^ - L 2 ^ _1 J) < a _<i • I det i?| _1 • 2 d • 2 l/3 l 

i=i 


Together we obtain 


/ h r (aBy) 1 dy = Y] / h r (aBy) 1 dy 

G « Gi 

= f p h r {aBy)~ l dy 

|/J|>dlog 2 a 2g “ 

< ^ 2 2r(d " |/3|) -a~ d -\ det B I" 1 • 2 d ■ 2 m 

|/3|>dlog 2 a 

oo 

< 2 2rd+d \detB\~ 1 -a~ d ^ 2 (1_2r)fc • |{/3 e Nq | |/3| = A;}| 

fc=[dlog 2 a] 


< 2 2rd+d \ det B\~ l ■ a~ d ^ 2 (1 " 2r)fc • (k + 1) 

!og2 a l 


d—1 


d—1 


dloff 9 a>l 

< ~~ 2 2rd+2d - 1 


= 2 2rd+d |det5|- 1 . a -^ 2 (1 - 2r)(fc+rdlog2 ^) • (fc + l+ rdlog 2 al) ,i " 1 

k=0 

oo 

< 2 2rd+d | det^l” 1 • a~ d ■ a (1 ~ 2r)d • ^ 2 (1 " 2r)fc • 2 d " 1 • (k + l) 4 *" 1 • |hilog 2 a] 

k =0 

oo / i \ d—l 

| det SI” 1 ■ a~ 2rd ■ 2 {1 ~ 2r)k (k + l)^ 1 ( 2d • ) 

fc=o ' °S / 

oo \ 

= I 2 2rd+3d ~ 2 d d ~ 1 \ det 5r 1 (log 2)^^ (2 1 “ 2r ) fc (jfe + l)*" 1 • a~ 2rd • (log a) 

k =0 / 

This is the desired estimate, since 2 1 ~ 2r < 1. 

Combining Theorem |T] and Lemma 0] yields: 


d—1 


□ 
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Theorem 2. Let B G M. dxd be a Frolov matrix and rgN. Then there is a constant 
c > 0 such that for every a > 2 1 ^ d and f G H r,mix (R d ) 


E|M 0 ,b(/)-/ d (/)| < c a~ rd ~ d ^ 2 (loga)^ 




The worst case error of M a> B for functions in H r,mix ([0. l] d ) is small, too. 

Theorem 3. Let B G M. dxd be a Frolov matrix and r G N. Then there is a constant 
c > 0 such that for every a > 2 1 ^ d and f G ZT r,mix ([0, l] d ) 

SU P I M a, B (f) ~ W) I ^ Ca ~ rd (log a) ^ ||/||ff7-,mix ([0i i ] d), 

LJ 

where the supremum is taken over all realizations Mf B of M a ^- 

Proof. The realizations Mf B of M a ^B take the form Q a nB.v for some u G [1, 2 1 ^ d ] d 
and v G [0, l] d . By Lemma [21 and Holder’s inequality, 


| QauB,v(f) — h{f)\ Z < I E | f(auBm) 

. mGZ d \{0} 


< I h r {auBm ) 1 

, mez d \{0} 


h r (auBm) ■ f(auBm ) 

K m£Z d \{ 0 } 


The first factor of this product is bounded above by a constant multiple of a~ 2rd ■ 
(loga) d_1 . This is proven similar to Lemma Q] 

Let N(/3) = {x G | |_2^' _1 J < \x 3 \ < 2^,1 < j < d} for j3 G Ng and 
Gf = {m G 7L d \ {0} | auBm G iV(/3)}. Then 7L d \ {0} is the disjoint union of 
all Gf over f3 G Nq. Again Gf is empty for |/3| < dlog 2 a. Otherwise, 


h r {auBm) > JJ (l + > Jj2 2r ^" 1) = 2 2r(l/3| - d) 

i=i j=i 


for m G Gf and hence h r (auBm) 1 < 2 2r ( d and 

2# 


iGfl < 


m G Z® \ {0} | |(Hm)jj < 


< 2“ !+l/3| a _ ‘ i + 1 < 2 d+1+l/3| a _<i , 
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since B is a Frolov matrix. This yields 


E hr(auBm)- 1 = E E h r {auBm)~ l < E 2 2r(cH/3 l } ■ a~ d ■ 2 d+1+ ^ 

mez d \{ 0} /36Ng mS G^ |/3|>dlog 2 a 

< ci • a~ 2rd ■ (loga) d_1 , 

like in Lemma [TJ 

We show that the second factor in the above inequality is bounded above by a 
constant multiple of ||/||^r,mix([ 0 yd)- This proves the theorem. For x G we have 


h r (x) ■ f(x) 


E 


c*e{0,...,r} d 


The function g a = D a f o ( auB) T has compact support in (aui?) T [0, l] d . Let 
M a = {/c G 7L d | k + [0, l] d fl (auB) T [ 0, l] rf ^ 0}. Then 


D a f(auBm ) 


D a f(y ) • e -^™Bm,y) dy 

' g a (x) ■ e~ 2ni{m ' x) da; 


det (auB) J Rd 
1 


det(aw.B) 






< 


\M a \ 


|det(ah5)| k£Mi 




L 2 (fc+[0,1]° 


Thus we obtain 


El h r (auBm ) • f(auBm) < El y] 

mGZ d \{0} mGZ d cte{0,...,r} d 

|M„| 


D a f(auBm) 


< 


| det (auB) \ 
\M a \ 

\det(auB)\“ 

\M a \ 

|det(awi?) | 


i E E E |<9<». e2 ’"'"‘'- , )L»( t+[ o..i 


<* feeM. 


ag{0,...,r} 

■ E E 

ae{0,...,r} d 

E n c ”/ii 


,2 _ |M a | 

' La (‘+M-) I det (auB)\‘ 


E NSc 


aG{0,...,r} d 


2 

L 2 (R d ) 


|M a 


| det (auB) \ 


H^ mi x([o,i] d ) • 


ctG{0,...,r} tl 

Since both \M a \ and |det(aw.B)| are of order a d , this yields the statement. □ 
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Error Bounds for H s (R d ) 

If, however, seN with s > d/2 and the integrand is from H s (JR d ) C C c (R d ), the 
following lemma holds. 

Lemma 5. For any Frolov matrix B G M, dxd and s G N with s > d/2 there is some 
c > 0 such that for each a > 0 and f G H s (R d ) 


'D a 


fix) 


da; < ca~ s+d/2 


H‘ 


Proof. Like in Lemma HJ we apply Holder’s inequality and get 


fix) dx^j = f(x) v s (x) 1/2 ■ v s (x) 1/2 dx 


'D a 


<( / v s (x) dx ) • \\f\\ Hsnad i < c 


'Da 


(Jd ( X + H^Hi) ' dx ) 


H°(R d ) > 


for some c > 0. Since ||x|| 2 > a for x G D a , the latter integral is bounded by 


-fxgR d : ||a;||2>aj- 


(l + ||x|||) s dx = / / (l + R 2 ) 6 • R d 1 dadi? 


« Sd -1 


= cr (Sd-i) / (l + i? 2 ) s ■ R d ~~ l dR < a (Sd-i) R~ 2s+d ~ l dR<c 


—2 s+d 


for some c > 0, since —2s + d — 1 < —1. □ 

In this case, combining Theorem [T] and Lemma [3] yields: 

Theorem 4. Let B G W lxd be a Frolov matrix, sgN and s > d/2. Then there is 
a constant c > 0 such that for every a > 0 and f G H s (R d ) 


E|-/„(/)!< ca— d/2 H/H^ . 

The Frolov property of B is important for Theorem 1 and the class LT r,mix ([0, l] d ), 
but we remark that B does not have to be a Frolov matrix to get this estimate on 
H s (/$L d ). As seen in the proof of Lemma 0 we do not need that the lattice points 
of aB1L d \{ 0} lie in D a but only that they he outside the ball {x G M. d j ||x ||2 < «}• 
For example, the identity matrix would do. But if B is a Frolov matrix, M a s 
works universally for H r,mix (R d ) and H s (R d ). Furthermore, the Frolov properties 
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of B prevent extremely large jumps of the number of nodes of M a ^B = QauB,v for 
small changes of a > 0 or u £ [1, 2 1 ^ d ] d . 

For functions from H s {[ 0, l] d ) the worst case error of M a ^ is also small. 

Theorem 5. Let B £ R dXcZ be a Frolov matrix and s £ N with s > d/2. Then 
there is a constant c > 0 such that for every a > 0 and f € LP([0, l] d ) 


su p| m £b(/)-W)| < ca s ||/||tfa([ 0) i]d), 

UJ 

where the supremum is taken over all realizations Mf B of M at B ■ 

Proof. The realizations Mf B of M a B take the form Q a uB,v f° r some u £ [1, 2 1 ^ d ] d 
and v £ [0, l] f/ . By Lemma [2] and Holder’s inequality, 

/ \ 2 


| QauB,v(f) - Id{f) | 2 < 


< I Y^ v s {auBm ) 

, mGZ d \{0} 


Y f(auBm) 
. mez d \{0} 


-1 


'Y. v s (auBm ) • f(auBm) 

Z d \{0} 


The first factor of this product is bounded above by a constant multiple of a 2s : 
Since 


v s (auBm) > \\auBm\\^f > a 2s ■ 11i3/77,||> a 2s ■ 11 B x | | 9 2a 


we have 


Yj v s (auBm ) 1 < a 2s ■ \ \B 1 

mGZ d \{ 0} 


Y W m 


mGZ d \{ 0} 


where this last series converges for 2s > d. 

We show that the second factor in the above inequality is bounded above by a 
constant multiple of ||/||^ s q 0 ip)- This proves the theorem. 

For any x £ M. d we have 


v a (x) 



Y. pvw 

IMIi<s 


The function g a = D a f o ( auB) T has compact support in (auB) T [ 0, l] d . Let 
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M a = [k G T> d | k + [0, l] d fl (omI?) T [0, l] d ^ 0}. Then 


D a f(auBm ) 


D a f(y) ■ 

' g a (x)-e- 2ni ^dx 


det (auB) J Rd 
1 


< 


det (auB) 
\M a \ 

| det (auB) | 2 




k£M a 






L 2 (fc+[0,1]° 


Thus we obtain 


v s (auBm ) • f(auBm ) < Z) |_D a /(cm.Bm) 

raEZ d \{0} mGZ d 

|M “' E E E |<9a,e“<"-->> ia(i+[ „, Ili; 12 


< 


mez d Nall, <skeMa 


\det(auB)\' 

\Mg\ II - _ _ 

|det(ahi?)| 2 Q " i2 ( fc+ M d ) |det(ah5)| 2 


\M a \ 


J2 ^ 


W 1<S 


| M a 


E ii fl “/n 


2 

L 2 (R d ) 


|M a 


|det(ahi?)| ^ J |det(ah£01"'' HS([ °’ 1]d) ' 

Since both |M a | and |det(aw.B)| are of order a d , this yields the statement. 


□ 


5 The Method M u g for Integration on H s ([0, l] <i ) 
and ff r - mix ([0, l] d ) 

We can transform the Monte Carlo method M a B from above such that it is still un¬ 
biased and its error satisfies the same upper bounds for the full spaces iJ T ' ,mix ([0, l] d ) 
and FP([0, l] d ), that M a)B satisfies for the subspaces i/ r,mix ([0, l] d ) and H s ([ 0, l] d ). 
This is done by a standard method, which is also used for deterministic quadrature 
rules for if r,mix ([0, l] d ), see 0 pp.359]. 

To that end let 0 : M —)• M be an infinitely differentiable function such that 
0|(-oo,0) = 0, 0|(i,oo) = 1 and 0|(o,i) : (0,1) -)• (0,1) is a diffeomorphism. For 
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example, we can choose 


h(x) 


e (2*-i)2-i if x e (0,1), 
0 else, 


if(x) 


5-00 h (t) dt 

IZ Kt) dt 


for i6l Like h also if is infinitely differentiable and obviously satisfies V’|(-oo, 0 ) = 
0 and ip |(i,oo) — 1- Since the derivative of if is strictly positive on (0,1), it is strictly 
increasing and a bijection of (0,1) with a smooth inverse function. 




with T(x) = (V>(ah), • ■ • ,if(x d )) T is 

-l(„\ _ \ \\T 


a diffeomorphism on (0, l) rf with inverse T 1 (x) = {if 1 (xi), ... ,if 1 (^d)) and 
|I?\l/(x)| detD^/(x) — Yli/j'(xi). 


i=l 

If Q is a linear quadrature formula for integration on the unit cube with nodes 
xJ J ) G [0, l] d and weights a 3 G R, where j = 1,... ,n, we dehne the transformed 
quadrature formula Q by choosing the nodes and weights to be 


and dj = aj ■ \D^f(x ^)|. 

Thus, Qs, v for v G R d and invertible S G R dxd takes the form 


Qs,v(f) 


n^si E f b ( S ~ T (™+<')))• h* ( S_T (“+0)1 

' ' ^ s-HI d 


for any function / : [0, l] rf —> R. Notice that |£)'L (S' T (m + v))| is zero, if 
S~ T {m + v) ^ [0, l] d . 

Algorithm. For a Frolov matrix B G R dxd and any a > 0 the transformed ran¬ 
domized Frolov quadrature formula M a> b is the method Q a uB,v with independent 
u and v, uniformly distributed in [1, 2 1 ^ d ] d and [0, l] d respectively. 


Lemma 6. The method M a B is well-defined and unbiased on L 1 ([0,1]^). 

Proof. Let / G L 1 ([0, l] d ). By the change of variables theorem ,/o = f ° V L • |DT 
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is also integrable on [0, l] rf and satisfies 


QauB,v(f ) = QauB,v(fo) and I d (f) = I d (fo)- 


Thus M“ B (f) = M“ B (f 0 ) for any realization Mf B of M a ^ B . This yields 

E (MaAfj) = E (MaAfo)) = Wo) = W) 


by Lemma [3l 


□ 


The following is our main result. It is important to recall that the number n 
of function evaluations in M aB is of the order a d , see Lemma [TJ 


Theorem 6. Let B £ M dxd be a Frolov matrix and r,s E N with s > d/2. Then 
there is a constant c > 0 such that for every a > 2 1//rf and f £ H r,mix ([ 0, l] d ) 


E 


M a , B {f ) 


sup 

L) 


KbV) 


W) 

W) 


< ca rd d/2 (loga)V H/|| ffr , mix([01]d) and 

< ca~ rd (loga)^ ||/||^, mix([01]d) 


and for every a > 0 and f £ H s ([ 0, l] d ) 


E 


M a , B {f) 


sup 

LJ 


MWf) 


W) 

W) 


< ca s d/2 ||/|| jKr - ([0>1] «* ) and 

— ca ll/llrr s ([o,i] ii ) > 


where the suprema are taken over all realizations Mf B ofM at B ■ 

Proof. Remember that Mf B {f) = Mf B (f 0 ) for / £ L 1 (M d ), / 0 = / ° E DT| 
and any realization Mf B of M aB . Since i/fx) = 0 for x ^ (0,1), we know that 
D a fo\d[o.i] d — 0 for each a £ {0,, r} d and hence f 0 £ H r,mix ([0, l] d ) C H r ’ mix (R d ) 
for each / £ Lf r,mix ([0, l] 01 ). 

This yields 


E 


M a , B (f)~W) =^\M a , B (fo)-Wo)\<c 


a 


—rd—d/2 


(loga) 2 • Wfo\ 


H r 


as well as 


sup 


m:,m)-w) 


= sup 




—rd 


(loga )^ 1 • ll/ol 


JJr, mix 
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if c > 0 is the maximum of the constants of Theorem [2] and Theorem [3j That 
proves the first statement, since there is a constant Co > 0 such that every function 
/ e H r ’ mix ([ 0, l] d ) satisfies ||/ 0 || fr r,m« (R d ) < c 0 \\f\\ H r^ m]dy 

This can be proven as follows. The partial derivatives of /o take the form 


D a f 0 (x) 


f}ll a lh 

dx • • • dx a d d 


d 

/('D) ■ 

l— 1 


Q lr-^ d 

£ D e Smx))-S a , e (x) 


Pi,—,Pd=0 


for a G {0,1,..., r} d , where S a ^{x) is a finite sum of finite products of terms 
^\xi) with i £ {1,... ,d},j G {1,..., rd + 1} and does not depend on /. ft is 
therefore continuous and bounded by some c a ^ > 0. Using the Cauchy inequality 


dimt; 

£’ 

2=1 


= |(U (1, ■ ■ ■ , 1) T )| 2 < Hulls • Hi 1 ) • • • > 1 ) T 112 = dimV ' 


TM 2 


dim v 


i —1 


for real vectors v, we get 




ir>"/ollV,< £ ||(z^/o»).s ojS | 


\Pi,...,p d =0 

«1 


L 2 ([0,l] d ) 




\Pi,...,p d =o 


a 1 ,...,a d 


C T(r+iy £ 

Pi,-,Pd=o 


= r 


= r 


I L 2 ([0,l] d ) 

ai,...,ad 

i)' £ 4,* / ^Kf^xyyix 


Pi,-,Pd=o 

ai,...,ad 


E 

a on J '3/ 


Pi,-,Pd=o 


f((0,l) d ) 

cti,...,ad 


|T/(T(T~ 1 (x))| 2 • |£T- 1 (a;)| dx 


<(r+l) d sup |D<f- 1 M| £ <^11^/11 


z£(0,l) d 

— c « ‘ ll/lln r > mix ([o,i] d ) > 


L 2 ([0,l] d ) 


Pl,...,Pd=0 


for some c a > 0 and 


||/o||//r,mix( R d) — || D /"o|| Z, 2 (M d ) — I ‘ 

a£{0,l,...,r} d \ a£{0,l,...,r} d 


m-.mixQo,l] d ) 
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The second statement is proven in the exact same manner. 


□ 


A translation of Theorem |6] by means of Lemma Q] shows that the algorithm 
M ay B indeed satisfies all the properties (PI) and (P2) and (P3). 


Summary. Let d,r,s £ N with s > d/2 and B £ R dxd be a Frolov matrix. Then 
there is a constant c > 0 (that may depend on B and r or s) such that for every 
n £ N there is some a n > 0 so that M anB uses at most n function values of any 
f £ L 1 ([0, l] d ) and satisfies 


E(M aniB (f) 


E 


ManMf) - W) 


sup 


<.,*(/)-W) 


E 


M^Mf) - W) 


sup 

UJ 


M^ B (f)-W) 


= W), 

< cn~ r ~ 1/2 (logn) V ||/|| Hr , mi x ([0il] d ) , if n>l, 

< cn" r (logn)V \\f\\ Hr ^ m]d) , if n>l, 

< cn~ s/d - 1/2 ||/|| Jf . ([0il]d) , and 

< cn ^ ll/lljysQQ^Jd) • 


Proof. Let c\ = (||5||i + I) 6 * > 1 be the constant of Lemma Q] and C 2 > 0 be 
the constant of Theorem |HJ For n > 4ci, we set a n = {n/{2c 1 )) l ^ d > 2 1 ^ d . By 
Lemma [U the Monte Carlo method M any B uses no more than 2ci • a 1 / = n function 
values of /. For n < 4ci we choose a n > 0 small enough such that the only node 
of M an B is zero. The method M anyB is thus unbiased for all n £ N and uses at 
most n function values. 

Theorem [6] yields for n > 4ci and thus a n > 2 1//rf that 


E 


E 


sup 

UJ 


sup 


(/) - W) 

< C 2 ( 

^ry d - d '\ 


C1>1 

< C 2 

d-^r (2cf) r+xl2 n- 


<c 2 ( 

(^rr~ d/2 wj 


= c 2 (2ci) s/<m/2 n - s / d - 1/2 


< C 2 ( 

(V 2 ci) 1/d )" rd ( log( 


log ( n / 2c 0 1/d 


lH s ([ 0 ,l] d ) 


d -1 
2 


Jf r >mix([0,l]rf) 
^r,mixQ0,l]d) ? 


d -1 
2 


H r, mix([0,l]d) 


d>l 


< c 2 d 2 (2ci) r n r (logn) 2 \\f \\ Hr ^ ixm]d) , 


k.,bU)-w) |< P 2 


^([o,i]h 


= c 2 (2ci) s/d n s/d ||/|| HS([0i i ]d ) • 
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This shows that the stated bounds hold for the maximum c of the constants 

C 2 <i 2 “ (2ci) r+1//2 , C 2 (2ci) S//d+1//2 and possibly larger constants that result from 

the cases n — 1,..., |_4ciJ. □ 
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